Unitary time-dependent superconvergent technique for pulse-driven quantum 

dynamics 
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I. INTRODUCTION 

The control of atomic and molecular dynamics by 
lasers has attracted considerable interest in the past 
decade. Time-dependent systems are traditionally stud- 
ied from a perturbative point of view with the Dyson 
expansion. The two limiting cases of a sudden and an 
adiabatic switching' of the perturbation have been ex- 
tensively studied In particular, this has resulted in 
the adiabatic theorem and the superadiabatic expansion 

aa 

We are interested in the case where the perturbation 
is switched on and off on a time scale which need not 
be arbitrarily small or large. It is well known that gen- 
erally, upon truncation, the Dyson series for the evolu- 
tion operator of a non-autonomous system is not unitary, 
giving rise to secular terms whose size grows with time. 
For periodic and quasi-periodic perturbations, a number 
of schemes have been proposed £|, notably one based 
on the Kolmogorov-Arnold-Moser (KAM) perturbation 
theory of classical mechanics 5j. Here we shall not be 
dealing with periodic or quasiperiodic systems but with 
time-dependent perturbations that are localized in time 
for which we develop, building on the KAM technique, a 
unitary superconvergent perturbation theory. 
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The KAM iterative method has been introduced in 
quantum mechanics by Belissard for periodic Hamil- 
tonians. It consists in generating at each step (with the 
help of a unitary transformation) a new reference or ef- 
fective Hamiltonian which collects higher order terms of 
the perturbation that commute with the reference Hamil- 
tonian constructed at the preceding step. At the first 
iteration the order of the perturbation is reduced from 
e to e 2 and, by considering the resulting Hamiltonian as 
a new starting point, the second transformation then re- 
duces the order of the perturbation from e 2 to e 4 . Hence, 
at the n-th iteration the size of the remaining perturba- 
tion is reduced from e 2 to e 2 . The quantum KAM 
technique has also been investigated for periodic Hamil- 
tonians by Combescure 0,0 an d more recently by Duclos 
and Stovfcek |9( . Quasiperiodic Hamiltonians have been 
considered by Blckher et al. in [Icj . 

All these authors have implemented the KAM algo- 
rithm in an extended Hilbert space constructed as the 
tensor product of the Hilbert space in which the origi- 
nal Hamiltonian is defined and the space of square inte- 
grable functions on the circle. This notion, introduced 
by Sambe in the periodic case and by Howland |T^ 
for more general time-dependent Hamiltonians, allows to 
construct a time-independent extended Hamiltonian (also 
called Floquet Hamiltonian in the periodic case) which 
is the starting point in the KAM algorithm. 

The KAM iterative procedure requires solving two 
commutator equations at each step. In |l3l IT^ Scherer 
has shown, adapting ideas from classical mechanics going 
back to Poincare, that these equations could be solved in 
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terms of time averages of some operators related to the 
perturbation. 

The generalisation of the KAM technique to time- 
dependent Hamiltonians has been worked out by Schcrer 
dm. It has been built in close analogy to classical 
mechanics and involves an extended phase space which 
includes time as a coordinate and the energy of exter- 
nal sources as its conjugate momentum, a notion closely 
related to that of [T^. On the other hand, the KAM 
algorithm proposed by Scherer is quite cumbersome to 
use and, in addition, is not guaranteed to yield a unitary 
evolution operator upon truncation. 

In this paper we present a KAM algorithm for non- 
autonomous Hamiltonians that is strictly unitary upon 
truncation at an arbitrary order. Moreover, for pulse- 
driven systems we construct explicitly the KAM trans- 
formations and study the convergence of the algorithm 
on a specific case. 

We start in Sec. Ill Al by recalling the KAM technique 
and the quantum averaging method for time- independent 
Hamiltonians. The notion of extended Hilbert space 
is presented in Sec. Ill HI at a purely formal level. In 
Sec. lII Cl we construct a unitary KAM algorithm for time- 
dependent systems in an extended Hilbert space. In 
Sec. Ill Dl the quantum averaging technique is extended 
to non-autonomous Hamiltonians in order to construct 
the KAM transformations directly in the original Hilbert 
space. Sec. IIII Al is devoted to perturbations that are 
switched on at some finite time in the past, for which 
we calculate explicitly the successive time averages in- 
volved in the KAM algorithm. We then focus on the 
case of pulse-driven two-levels systems, for which we re- 
sum exactly in Sec. IIII Bl the infinite series of commu- 
tators yielding the remaining perturbations at a given 
step of the iterative procedure. Finally, in Sec. IIII CI the 
method is applied to a two-level system interacting with 
a sine-squared pulse, taking the ratio of the characteristic 
duration of the pulse and the characteristic time of the 
free evolution as the small parameter e. We show the re- 
markable result that the KAM algorithm converges for all 
values of the parameter e, even larger than unity, allowing 
to go from the sudden regime to the opposite adiabatic 
regime. The conclusions are given in Sec. II VI while some 
details of the calculations are reported in Appendices A 
and B. 



II. UNITARY SUPERCONVERGENT 
TIME-DEPENDENT PERTURBATION THEORY 

A. KAM algorithm for autonomous Hamiltonians 

In this section, we present the KAM technique for 
a time-independent Hamiltonian following the formula- 
tion of and using the averaging method of Q. Let 
K\ = Ko + eVi where Kq is a reference Hamiltonian de- 
fined on a Hilbert space H and eVi a bounded self-adjoint 
perturbation with small parameter e. As will become 



clear below, the subscript n indicates that an operator 
A n is involved in the n-th iteration. On the other hand, 
the upperscript e stands for effective and indicates that 
an operator constructed at the n-th step will be taken 
as the new reference at the next step. Throughout the 
paper the leading order in e will appear explicitly in front 
of the operators which are thus themselves of order e° but 
may still depend on e although we shall not indicate it 
explicitly. We look for the generator W\ of a unitary 
transform T\ = e eWl such that 

T\K X T X = T\ (Ko + eV 1 )T t = K\ + e 2 V 2 = K 2 , (1) 

with [K^, K ] = 0. Writing K\ = K + eDi, the unknown 
W± and Di are solutions of the following commutator 
equations: 



[#o,A] = 0, 

[K ,W 1 } + V 1 =D 1 



(2a) 
(2b) 



The remainder e 2 V 2 contains all the terms of Eq. 
which are not of order e° (which disappear trivially) or of 
order e (which disappear identically because of Eq. I[2B|| ). 
It reads 



e 2 V 2 



2 2 

\[W 1 ,V 1 ]- e -[W 1 ,D 1 ] 

3 3 

+ j[W 1 ,[W 1 ,V 1 ]] + j[W 1 ,[W 1 ,D 1 ]] + ..., (3) 

or, writing the series of commutator in a compact form 
that we shall use later, 

Jy a = E 7 fc + lV {k^iWuV^ + ad^WuDt)], 

k=l ^ 

(4) 



where 



D 



k = 



ad(A,B) = i \ A&d k-i iA)B) \ fc >j 



(5) 



The solutions to Eqs. © can be written in terms of av- 
erages: 

l-T 



lira — 



UK, 



°V 1 e ltKo dt 



A'n 



(6a) 



[ dt [ dt'e- lt ' Ko (Vi 



-m 



~K 



(6b) 



This is readily checked upon substitution, noting that 
e -itK Q j g -j-j^g p r0 p a g a tor of the reference Hamiltonian Kq 
in units such that h=l. The shorthand notations on the 
right handside of Eqs. © can be viewed as well defined 
linear transformations of the operator V\. 

The process can be iterated, transforming now the op- 
erator K 2 defined by Eq. JQl with T 2 = e e W2 and con- 
sidering K1 as the new reference operator: 

T\ K 2 T 2 = Tl{Kl + e 2 V 2 )T 2 = K c 2 + e 4 V 3 = K 3 , (7) 
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with [if |, Kf] = 0. Notice that the new perturbation 
is not of order e 3 as would be the case in a standard 
perturbation theory, but of order e 4 . Similarly, after n 
iterations we obtain 

T 7 \K n T n = if I + e 2 " V n+1 = K n+1 , (8) 

with [if^if^-J = 0. The new reference Hamiltonian 
if ° is written 

n 

Kl = K e n _ x + e 2 ^D n = if + Y, ^ D o ■ ( 9 ) 

i=i 

The generator W n of the unitary transformation T n = 

e e Wn and the operator D n , solving equations analo- 
gous to Eqs. are calculated as 

D„ = Vf 1 - 1 , (10a) 

W n = -iVn*^ 1 . (10b) 

The remainder being of order e 2 ™ , the KAM algorithm 
is called superconvergent. However, we emphasize that 
a proof of convergence is to be established in each case. 
We also note that in the absence of resonances, e need 
not be small for this algorithm to converge ^3] • 

B. Extended Hilbert space 



H(t + s) is considered as an intermediate step). The lift 
of Eqs. (| 1 1 11 on the extended Hilbert space K reads 

^U n {t, to; a) = H(t + s)U n (t, t ;s) , (12a) 
Uu{t ,h\s) = tj J ®±u- (12b) 

where 1l is the identity operator on L. Notice that 

U H (t + s',t + s';s~s') =U n (t,t ;s) Vs' e R. (13) 

Finally, an extended Hamiltonian is defined on K as 
the time-independent self-adjoint operator K, = Ti.(s) — 
i-^ <g> 1h- Its associated unitary evolution operator reads 
Mic(t,to) = e - ^* - * ^ and is related to the solution of 
Eqs. (|12fl by the following equation which is easily de- 
rived: 

W K (Mo) = T- t UH(t,t ;s)T to , (14) 

where the translation operator T* acts on functions 
£(s) G L according to Tt£(s) = £(s + t) and can be ex- 
pressed as T* = e t- B*. 



C. KAM algorithm in the extended Hilbert space 
for non-autonomous Hamiltonians 



Given a time-dependent Hamiltonian H(t) acting on 
a Hilbert space H, we recall here how the notion of ex- 
tended Hilbert space of [12j allows to construct a time- 
independent operator on that space. Let Unit^to) de- 
note the evolution operator associated to H(t), so that 
the Schrodinger equation and the initial condition read 

i^-U H (t,t ) = H(t)U H (t,t ) , U H (t ,t ) = t m , (11) 
at 

where In is the identity operator on H. We introduce 
a parameter s £ R which plays the role of an arbi- 
trary reference time, and let Uu(t,to; s) be the solu- 
tion of Eqs. (|llfl now with H(t + s). In H, the opera- 
tor Unit, to; s) depends parametrically on s. Notice that 
U H (t,to;0) = U H (t,to). 

An extended Hilbert space K where s is now an ad- 
ditional coordinate can be defined as the tensor product 
of L and HI where L = £ 2 (K) is the space of square in- 
tegrable functions on the real line: K = L <g> H. The 
family of operators Uii{t,to] s) acting on H is lifted to 
the operator Unit, to; s) defined on K by considering the 
full dependence on s as a multiplication operator on L. 
Similarly, the family of operators H(t + s) on Et is lifted 
to the operator H(t + s) on K. We shall denote operators 
acting on the extended Hilbert space IK by uppercase cal- 
ligraphic letters and shall refer to U-}{(t, to; s) and H(t+s) 
as the lifts of Uu(t,to) and H(t) respectively (with the 
understanding that the family of operators Un(t, to; s) or 



We consider eVi(i) as a bounded time-dependent per- 
turbation of the time-dependent reference Hamiltonian 
H (t) defined on the Hilbert space H and whose propa- 
gator UH (t,to) is known. Our aim is to obtain a KAM 
expansion for the evolution operator (t, to) of the full 
Hamiltonian H\(t) = H (t) + eVi(t). We first consider 
the extended Hilbert space K = L ® H, and the lifts 
^■Hi(Mo;s) and TL\{t + s) of the operators fTjj^i, io) 
and iii(t) as defined in Sec. Ill HI Similarly, Un (t,to;s), 
Ti.o(t + s) and Vi(t + s) denote the lifts on K of UH a (t,to), 
H (t) and V\{t). We then define the associated extended 
Hamiltonian on K: 



JCx=Hx{s) -i— ® 1 H , (15a) 
os 

= H (s) - O 1 H + eVi(s) = /Co + eVi(s), (15b) 
os 

which is of the form considered in Sec. Ill Al Hence, 
we can now apply the KAM technique in the extended 
Hilbert space to obtain a KAM expansion for the time- 
independent operator Kx- 

At the n-th iteration of the algorithm, the operator 
K n = K, c n _ l + e 2 V n (s) is transformed by the unitary 
operator T n (s) according to 

T n \s)K n T n {s) =K c n + e 2 "V n+1 (s) , (16) 
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with The remainder e 2 V„+i (s) is given by an expression anal- 

ogous to Eq. 0: 

*2 = K(^)-^®1h, (17a) 



KM = K-iW + ^Vn "-'(a), (17b) 
T„(s) = exp^-ie 2 "" 1 ^" -1 ^)^ . (17c) 



e 2 "V„+i( S ) = £ 



fc=i 



jfe e (fe+l)2" 

1)! 



fcad K V„ " _1 (s),V„( S ) + ad K V n " (s),V„ "^(s) 



(18) 



The operators Vn*'™ 1 (s) and V„ ™ 1 (s), defined by Eqs. ©, can be expressed in terms of the operator Un e _ (t, to; s) 
using Eq. (fHfjl for the propagator of /Q : 



V n K - 1 (s)= lim i 

i — >oo J 



dtW K »_ 1 (0, -t; s)V n (s - t)Z^._ i (0, -t; s) = V„ ""'(a), 



(19a) 



^"'(s) = lim- / di'/ diW W c (0,-t; S )\v n ( S ~t)^V n n - 1 {s~t)]ul ic (0, -t; s) = («). (19b) 



JO 



On the other hand, from Eqs. I|14|) and l|17() and the fact that K, c n commutes with all the operators K.% constructed at 
the preceding iterations one deduces that 



W«e (t,to; s) = Un (t, t ; s) Si(*,*o; s + t Q ) . . .S n (t 7 t ; s + to), 
= Si (t, t Q ;s + t) ... S n (t, t ;s + t) U Ho (*> t ;s), 

I 



(20a) 
(20b) 



where S p (t,to; s) denotes the following operator on K: 
S p (t, t ; s) = exp (-»(* - t^e 2 "" 1 ^- 1 (*)) . (21) 

The detailled derivation of Eqs. (|19|l and lj20|l is provided 

in appendix 1X1 It follows that V„ nl (s) and V„ " (s) 
are calculated from W^ (t, ioj s) as well as the operators 

V p p-1 (s) and V p P ~ (s) constructed at the preceding 
iterations. Hence, the operators T n (s), K? n and V n +i{s) 
entering Eq. (|16f> are now entirely determined. 

The extended Hamiltonian JC\ we started with can be 
expressed in terms of /C° by repeated use of Eq. I|16|) : 

K x = Tiia) . ..T n (8)K? n t[(8) . ..T?(s) + 0(e 2 "). (22) 

The propagator WjCi(* } *o) = eT^-M^i a ll ows then to 
construct the operator Un 1 (t,to; s) on K from Eq. (|14|l . 
Taking also Eq. (|20a|l into account yields 

U-Ht (t, t ; s) = Tx (t + s) . . . T n (t + s) U Ho (t,t Q ;s) 
Si(t,t ; s + to) . . .S n (t,t ; s + t ) 
Zl(t + s)...Tl(to + s)+0( ( ? n ). (23) 

It is now possible to return to the original Hilbert space 
H, by considering the dependence on the s— variable of 
each of the operators entering Eq. (|23|) . which defines a 



multiplication operator on L, as a parametric dependence 
on time in H, and subsequently setting s = 0. Hence, in 
agreement with our notations, Eq. ill'-'il) is the lift on IK 
of the following expression for the propagator Uuxit^o) 
on H: 

U Hl {t, t ) = Ti(t) . . . T n (t) U Ho (t, t )Sx(t, t ; t Q )... 

S n (t, t ; t ) Tt(to) . . . T}(t ) + 0(e 2 "), (24) 

where S p (t, to; to) and T p (t) acting on H are obtained, as 
just described, from their lift S p (t, to; s + to) and T p (t + s) 
constructed on K. 

In practice, however, we shall find it simpler to con- 
struct S p (t,to;to) and T p (t) directly in H. As we show 
below, this can be achieved by considering Eqs. I|17b|) - 
(|21|l which have well defined meaning on K as the lift of 
equations for corresponding operators defined on H. In 
particular, Eq. I|19a|l with n = 1 is the lift of a similar 

equation defining the operator V\ H ° {t) on H in terms of 
U Ho (t,t ) andy x (t). 

D. KAM expansion in the original Hilbert space 
for non-autonomous evolution operators 

In this section, we shall construct the propagator 
Uii 1 (t,to) from Eq. (|24|) through an iterative procedure 
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entirely denned in the Hilbert space H of the Hamilto- 
nian Hi(t) = H Q (t) + eVi(t), i.e. we shall not have to 
define operators in an extended Hilbert space. The oper- 
ators UH (t, to) and V\(t) allow to construct the operator 

Vi ° (t) on H according to Eq. (|26a|l given below, where 

we set p = 1. Subsequently, the operator V\ ° (t) can 
be obtained from these operators using Eq. (|26b(l with 
p = 1. Hence, by Eqs. (J2SJ, the operators Si(t,t ;t ) 
and T±(t) are determined. 

On the other hand, Eq. l|26c|l with p = 1 enables us 
to derive the operator V2(t) on EI from the operators 

V\ H ° (f) and Vi °(t) we have just constructed. Similarly 



Eq. (|27|l yields the propagator Unfit, to)- It follows that 

the operators V^ 1 ft) and V2 1 (i) can be obtained from 
Eqs. (|26a() and (|26b|l now with p = 2. 

For p > 1 we have the following operators on H: 

$,(<, to; ti) = exp (-»(< - io)e 2 ^ 1 ^- 1 (ti)) , (25a) 
T p (t) = exp (-ie^'Vp"'- 1 (t)j , (25b) 

where 



V P B '- I (jt) = lim 



1 



1 



du tfflj^ U )^,(u)t^j_ i (t, u), 



t-T 
T 



% p -' (t) = lim - / dt' du U H; (t, u) \v p (u) - V/**- 1 (u)] (t, it), 

00 -ft (fc-l)2 p_1 r / H c \ / H c h<= 

V p+1 (t)=J2 ru 1 im 1 fcadfc K P_1 (*),^W) +ad fe (^ "-'(t),^*- 1 ^) 



fe=i 



(fc + 1)! 



(26a) 
(26b) 
(26c) 



and 

U H} (t,t ) = f/ Ho (Mo) Si (Mo; to) • ■ ■ S p (Mo;*o). (27) 

Note that = F«_ x (t) + e 2 ^ 1 V^" 1 (t). By the it- 

erative procedure described here and which rests solely 
on Eqs. the operators S p (t,t ;t ) and T p (t) 

are constructed entirely in the Hilbert space EL The 
propagator [/^(Mo) of the non- autonomous Hamilto- 
nian H x (t) = H (t) + eV\{t) is then obtained by Eq. 
up to a desired order in e. 



III. PULSE-DRIVEN SYSTEMS 

A. General case 

In this section, we consider the physically relevant 
case of time-dependent perturbations eVq(t) which are 
switched on at a given finite time ti. To allow for some 
flexibility in the choice of the reference operators we shall 
consider the slightly more general case of perturbations 
eVi (t) which before ti are constant in time and commute 
with the reference propagator Z7# (Mo) on HL 

V 1 (t) = V 1 (t i ) V*<ti, (28a) 
[V 1 (t i ),U Ho (t',t)] = Vt,t'<ti. (28b) 

After ti the time-dependence of Vi(t) is supposed to be 
uniformly bounded in time but ortherwise arbitrary. In 
particular, it need not be turned on or off infinitely slowly 



or rapidly, and need not be constant or periodic in the 
meantime. For this class of perturbations, that we refer 
to as pulsed perturbations, the limits in Eqs. I|26|l can be 
calculated as we show in appendix IBI 

On the one hand, from Eq. (|26a|l we obtain 

V[ H °(t) = U Ha {t,ti)Vi(ti)U^ Ho {t,ti), (29a) 
Vp H Z-i(t)=0 V;p>l. (29b) 
Hence, Eq. (|25a|l yields 



S x (t, t Q ; t ) = U Ho (to, t-Je-W-^WWuh (to, U), (30a) 
S p (t,t ;t ) = t M V;p>l, (30b) 

which by Eq. (|27|l implies 

U H «(t,to) = C/ ffo (Mi)e- j( *- to)eV ' l(t ' ) C/; o (to,ti)- (31) 

Note that ff^_i(i) = Hf(t) for all p > 1. 
On the other hand, Eq. I|26b|) results in 

Vf°(t) = J duU Ho (t,u) [Fi(w) -TT H °(u)] C4 (*,«). 

(32a) 

t 

dul/fl;(t,u)^(u)E/& ! (t,u) Vp > 1, 

(32b) 

where Vp(t) is given by Eq. (|26c|) . 



v p Hp - x (t) 
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For pulsed perturbations it follows that Eq. I|24|) for 
the propagator Uh x (t, t ) up to correction terms of order 
e 2 becomes 



U Hl (t, t ) =T 1 (t)... T n (t) U Ho (t, f i ) e -i(t-*o)eV 1 (t 1 ) 
uUt ,t i )Tl(t Q )...Tl(t Q ) + O(e 2n ), 



(33) 



where T p (t) defined by Eq. (|25b|) is obtained through the 
simple integral given in Eqs. 13 2() . 



B. Exact resummation of the remainder e 2P V p +i(t) 
for pulse-driven two-level systems 

For two-level systems some of the formulas that we 
have constructed in the preceding sections can be writ- 
ten in an explicit simple form. In particular, we shall 
calculate exactly the remainders of the KAM iterations. 
The partition of an Hamiltonian Hi(t) = Ho(t) + eVi(t) 
on HI = C 2 can always be choosen such that 



Vi(t) =^v k {t)a k , 



(34) 



fe=i 



where v k (t) G £2^) are real functions on R and a k the 
Pauli matrices: 



0"! 



1 

1 



(7 2 



-i 

1 



0"3 



1 

-1 



(35) 



The perturbation eV\ (t) is switched on at a finite time i; 

so that Eqs. l(2^|) reduce to V p H ° (t) — for all p > 1, and 
Eqs. (|22J) to 

V p H \t) = f duU Ho (t,u)V p {u)U^ (t,u) Vp > 1. (36) 



Furthermore, the infinite series of Eq. (|26c() for V p +i(t) 
with p > 1 reads 



v p+1 (t)=j2 



fci fc e (fc-l)2» 



( fc + 1 ) ! 



ad' 



(y p H °(t),V p (t)). (37) 



Let B p (t) be the unitary matrix which diagonalizes 

V p (t). As the matrix V\{t) is traceless, it is straight- 
forward to show by induction with the help of Eqs. i|36l) 

and H37fl that V p ° (t) is traceless for all p > 1. Hence 



Bl(t)V p (t)B p (t) = -r p (t)a 3 , 



(38) 



where 



r p (t) = [-detV p Ho (t)}^ 2 , (39) 
is real. The following identity holds for all k > 1: 

ad fe (y p H ° , V p ^ = B p ad k (-r p a 3 ,BlV p B p ) B^ (40) 
and for any Hermitian matrix M in C 2 one has 

ad 3 (cr 3 ,Af) = 4ad 1 (cr 3 ,M). (41) 

Combining Eqs. (pulf and (f^Tf) yields 

(V P H ° , V p ) = (-2r p ) fe - W (V P H ° , V p ) , (42) 



ad* 



where £ is 1 if k is odd, and 2 if k is even. The series of 
Eq. H37|) for Vp+i(f) can then be cast into the form 



v P +i(t) = [vp ffo (;),v p (i)l +e 2P " 7p (i) k ff °(t), k ffo (t),vi(t) 



where in agreement with our notations the following quantity are of order e°: 



7p(*) = 



.cos[2e 2P r p (f)] - 1 + 2e 2P r p (t) sin[2e 2P r p (i)] 
* [2e 2p -V p (t)] 2 : 

2e 2 ^\ p {t) cos[2e 2P "Y p (£)] - sin[2e 2P ~V p (t)] 
[2e^r p {t)Y 

I 



(43) 

(44a) 
(44b) 



The remainder e 2P V p +i(t) is therefore well-defined for all 
values of the parameter e and all times t. We then ob- 
tain from Eq. I|33fl the propagator Uu^i, to) up to an 



arbitrary order in e, 



U Hl (t,t )=U$(t,t ) + O(e 2 '") 



2" - 



(45) 
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where 

U$(t,t Q ) = Ti(t) . . .T n (t)U Ho (t,t )Tl(t ) . ..T^to). 



(46) 



Moreover, using Eqs. I|25b|) and (|38[1 one deduces that 



T p (i)=cos[e 2 r„(t)]l 



c 2 



r„(t) 



V P (t), 



-H 



(47) 

where V p (t) is given by Eq. (|36l) , and owing to Eq. 143(1 
can be calculated for arbitrary p without having to re- 
sort to an infinite series. Note that the form of Eqs. I|44|) 
suggests that Eqs. ()43[) and (|47|l remain valid for finite 
values of the parameter e, possibly larger than unity. We 
shall see below that this is indeed what is observed nu- 
merically. 



C. Numerical implementation for a two-level 

system 

Here we investigate the convergence of the KAM tech- 
nique with the number of iterations as well as its domain 
of validity for a specific two-level model perturbed by 
a pulsed interaction. The algorithm is implemented nu- 
merically for a system described by the time-independent 
Hamiltonian loo~3 and which interacts through a\ with 
a sine-squared pulse of characteristic duration r. This 
pulse shape is commonly used in the literature because 
of its bounded support and continuous first derivative at 
the boundaries. Defining the characteristic duration r as 
twice the full width at half maximum fixes the total dura- 
tion of a cycle to r and yields the following dimcnsionlcss 
pulse shape between the dimensionless time t% = and 
t f = 1: 



Q(t) 



2Asin 2 (7rf) 




< t < 1, 
elsewhere. 



(48) 



Note that the peak amplitude is twice the pulse area 
A = A(t{) where A(t) = f. Q(u)du, and that it can 
be fixed independently of the parameter e = ujt that 
we shall take here as the small parameter. This allows, 
in particular, to treat large non-perturbative areas for 
short pulse durations, which corresponds to the exper- 
imental conditions used to generate short laser pulses. 
The Schrodinger equation reads 

i^U{t, t ) = [n(t)o-i + e<j 3 ] U(t, t ), (49) 
with U(to,to) — Ic 2 - For e = its solution is 

£/ (0) (Mo) = e -^«-^(*°)K (50) 

As a first step, we write Eq. (|49|l in the interaction 
representation with the help of the unitary operator 
(7 (0) (Mi): 

i^U Hl (t,to) = HtitpH^to), (51) 



where 

fli(i) = eU^fatfiaaUWfati), (52a) 
U Hl (t,t ) = D'( ) t (t,ti)U-(t,to)^(to,ti). (52b) 

Note that .ffi(ti) = £03 so that if we substract £173 from 
Eq. (|52a|l we obtain an operator that vanishes for t = t 1 . 
It is interesting to identify the latter as the perturba- 
tion eVi(t) in order to get rid of the average V\ Ha {t) [cf 
Eq. (|29al) ] as was done in Sec. IIII Bl 

VL(t) = U<f>'> 1 (t,t l )<r3U ( ®(t,t l )-<r3, 

= (cos[2A(t)] - 1) 03 + sm[2A(t)]a 2 - (53) 



The reference operators are thus 



H 

U Ho (t, t ) 



-i(t— t )ecr 3 



(54a) 
(54b) 



The KAM algorithm is now applied to Eq. JHH. Af- 
ter n iterations this results in Eq. I|45|l . The propagator 
U(t,to) is then obtained from Eq. I|52b|l : 



U(t,to) = U {n \t,to) + 0(e 2 "), 



(55) 



where 



U^(t,t ) ee UM(t,W$(t,t )UW 1 (to,td, (56) 

in which uj$ (t, to) is given by Eq. |@5J. 

Let | + 1) and | — 1) denote the first and second column 
of g 3 respectively. For given e and A, the wave function 
at the time the pulse is switched on is tp(ti) = \ — 1). At 
the end of the pulse, the error between the wave func- 
tion ij)(tf) = U(tf,ti)ip(ti) computed by solving numer- 
ically the Schrodinger equation and the wave function 
-0(™)(£ f ) = C/C") (t f , ii)V(ii) obtained after n KAM iter- 
ations (where the integration in Eq. I|3()|) is performed 
numerically) is defined as 



E |toh%)>-v7lV> (n) (*f)> 



.7J = ±1 



1/2 



(57) 



Figure^ displays the accuracy of the KAM algorithm 
for different number of iterations as a function of e in 
the case where the pulse area A = | , which corresponds 
to a peak amplitude equal to tt. One can see that this 
procedure reproduces the numerical results with great ac- 
curacy for any value of e provided a sufficient (yet small) 
number of iterations is used. For the range of e shown 
on Fig. Q the fifth iteration is indistinguishable from the 
numerical result. For e = 0, the peak amplitude consid- 
ered here leads to a complete population transfer from 
the lower to the upper state (the so-called "7r— pulse" 
transfer) . Figure ^3 shows that this transfer decreases 
for larger e until becoming negligible beyond e w 5, a 
feature which characterizes the adiabatic regime. It is 
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For A = 5 and different number of KAM iterations 



FIG. 1 

n, (a) natural logarithm of the error Za n defined in Eq. l|57). 
and (b) transition probabilities |(+1|t/j)| 2 computed numeri- 
cally (thick line) and |{+l|V> (n) )| 2 calculated with the KAM 
algorithm (thin lines) as a function of e, at the end of the 
pulse. 



striking that the adiabatic regime can be reached with 
great accuracy from the third iteration on. 

The accuracy of the KAM algorithm is plotted as a 
function of the number of iterations in Figure [21 As 
suggested by the order e 2 of the remainder, the error 
decreases faster than exponentially. 

Figure [3] displays the accuracy of the KAM algorithm 
as a function of the pulse area A. As expected by in- 
spection of the Schrodinger equation l|49|l . Fig. \Qp shows 
that for larger pulse area the pulse is effectively more sud- 
den, since the transition probability can reach maximum 
values closer to 1. The KAM algorithm accuracy is con- 
sequently globally better, except for pulse area smaller 
than 7r/2, as seen on Fig. 



IV. CONCLUSION 

We have derived a unitary superconvergent algorithm, 
based on the KAM technique, that allows to treat time- 



dependent perturbations that are localized in time. In 
the physically relevant case of perturbations that are 
switched on at some finite time in the past, we have 
shown that the computation of the KAM transforma- 
tions can be greatly simplified. The remarkable efficiency 
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FIG. 2: Natural logarithm of the error A n as a function of the 
number n of KAM iterations for A = -| and e = 0.5 (crosses), 
1 (stars) and 2 (circles). 



of the method has been shown for a pulse-driven two- 
level system, for which we obtain convergence all the way 
from the sudden regime to the opposite adiabatic regime. 
We anticipate interesting applications of this method in 
the context of alignment and orientation of molecules by 
pulsed laser fields. 
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APPENDIX A: KAM ALGORITHM IN THE 
EXTENDED HILBERT SPACE 

1. V K {s) and V K (s) 

For time-dependent problems, the KAM algorithm in- 
volves calculating the following transforms of operators 
V(s) with respect to the propagator of JC = H(s) — i-i- ® 



1h on the extended Hilbert space 



Li 



V K (s)= lim i / dt 

v*00 



- i " c V(s)e i *' c . 



lim — 



t' 



T 

dt' I dt e 
o Jo 



-UK 



(V(s) 



V K ( S 



(Ala) 



UK 



(Alb) 



Hence, one has to consider operators on K of the form 
B(s,t) 



-UK 



V - V K for V K . Using Eq ifTi)). B(s, t) becomes 

B(«,t) = T_tZ/ w (i,0;s)^(«)^(t,Q;B)T t 

= W w (t, 0; s - i)A(s - 0; s - t) 

= U n (0 : -t; s)A(s - t)lth(P, -t; s) , (A2) 

where we used the definition of the translation operator 
and Eq. (|13|l to obtain the second and third equalities, 
respectively. 

At the n-th iteration of the KAM algorithm, Eqs. lfT?|) 
imply taking IC = fCt^ and V = V n in Eqs. (|A l|l . hence 
U H = Un° n _ 1 in Eq. i|Al)l. 
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2. Un-Jt,t ;s) 

We show here that the operator Uu^ (t, to; s) can be cal- 
culated according to Eq. (|20a|l . or equivalently Eq. (|20b|) . 
in terms of UH (t,to) and the operators Sk(t,to', s) with 
1 < k < n defined by Eq. J2J . Indeed, applying Eq. itTHl 
to the propagator of /C® defined by Eqs. ifTTjl yields 

Uni(t,t ;s) = TtU^i (t, t )T- to 

= T t Uici 1 (t,to)S n (t,to;s)T-t 

= T t ^e_ i (t,to)T_ to 5 n (t,t ; s + t ) 

= U H ^_ t (t, t ; s)S n (t, t ; 8 + to), (A3) 

which proves Eq. J3UaJ. Note that by construction 

[K-n-iiVn ™ _1 ( s )] = 0> which is crucial for writing the 
second equality. This commutation relation also allows 
to permute i and S n on the second line of Eqs. 1A3L 
which then results in Eq. (|20bl) . 

APPENDIX B: AVERAGING FOR 
PULSE-DRIVEN SYSTEMS 

In this appendix, we consider the case of a time- 
dependent operator V(t) which, before some finite time 
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ti, is constant in time and commutes with the propagator 
of an Hamiltonian H(t) on H: 

V(t) = V(U) Vt < t u (Bla) 
[V(ti),U H (t,t )] = Vt,t <*i- (Bib) 

After ti the dependence on time of V(t) is arbitrary pro- 
vided it is uniformly bounded. We show that the op- 
erators V (t) and V H (t), defined by Eqs. I|2tj[l . can be 
calculated as 



where 



(B2a) 



V H (t) = J AuU H {t,u) \ V{u)-V (u)J Ujj(t,u).(B2b) 



We first prove Eq. l|B2a|) . rewriting Eq. (|26a[l as 

V H (t) = lim - / duA(t,u), (B3) 

T^oo J J t _ T 

where A(t,u) = Un{t, u )V{ u )Uh (t, u). The propagator 
UH(t,u) can be decomposed into u). If 

u < ti then Uh(U, u) and V(u) satisfy Eqs. IjBlfl implying 



A(t, u) = I7jEr(t, ti) < ti. 



(B4) 



For f < t; the integration domain in Eq. (|B3jl is such that 
A(t, u) is given by Eq. (|B4fl . resulting thus in Eq. {B2aJ. 
Notice that this latter reduces to 



-H 



V (t) = V(ti) Vt < ti. 



(B5) 



On the other hand, when t > ti there is a remaining 
integral over «£ [ti, t] which is bounded and independent 
of T, vanishing thus in the limit. The result given in 

Eq. lfB2ajl follows. 

We now turn to the proof of Eq. (|B2b|l and write 
Eq. l|26bTl as 

V H {t) = lim i / dt' [ duB(t,u), (B6) 

T^co T J Q J t _ t , 



B(t,u) = U H (t,u) V{u)~V (u) U]j{t,u). (B7) 



The case t < ti is directly proven since B(t, u), which is 
also the integrand in Eq. !B2b|) . vanishes identically by 
Eq. JB5J: 



(B8) 



V H (t) = vt < u. 



For t > ti splitting the domains of integration in Eq. I|B6 
yields 



V H (t) = lim^ i U dt'J^duB(t, u) 

+ [ dt' [ duB(t,u)+ f dt' [ duB(t,u)\. (B9) 
Jt-ti Jt-t' Jo Jti J 

The first double integral being bounded and indepen- 
dent of T does not contribute in the limit whereas the 
second one vanishes because B(t,u) = if u < ti. The 
result of Eq. lB2b|l comes from the last double integral 
of Eq. El . 



Finally, we show that if Eqs. I|B1|I are satisfied with 
V = Vx and H = H , then Eqs. (B2J hold with V = 
V n and H = H°_ 1 for any n > 1. The case n = 1 
follows directly. We prove the case n > 1 by induction, 
assuming Eqs. I|B1|) are verified for n — 1 . The operator 
V n (t) is obtained by Eq. (|26cll as a sum of terms involving 
-ff°_ 2 

V n -i (t), which by Eq. {BBJ) is zero for t < ti. Hence 
V n {t) = for t < U so that Eqs. (JbTJ hold for n > 
1, which concludes the proof. Notice that for n > 1 
Eq. fB2aT yields Vn" 1 ' 1 {s) = 0. 



